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Discovery of a galactic wind in the central region of M100 
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ABSTRACT 

We report the discovery of a galactic wind in the central region of the galaxy M100. 
This result is based on a careful 2D spectroscopic study performed on observations 
made with the fibre system INTEGRAL on the WHT. The primary evidence of the 
wind is the presence of blueshifted interstellar NaD absorption lines. The velocity field 
of the absorbers show a clear rotation pattern but globally blueshifted (<~ -115 km/s) 
with respect to the systemic velocity of the galaxy. The emission lines also present 
a blueward component arising from the ionized gas phase of the galactic wind. The 
velocity field of the ionized gas wind component shows no evidences of rotation but 
exhibits a pattern that can be interpreted in terms of the projection of an outflowing 
cone or shell. The wind component has [Nil] /Ha ratios of about 1.8, typical of shock 
ionization. The ionized component of the wind can be identified with an expanding 
shell of shocked gas, and the neutral component with disk gas entrained in the wind 
at the interface of the expanding shell with the galactic ISM. The galactic wind seems 
to be driven uniquely by the nuclear starburst. Our analysis indicates that a non 
negligible fraction of the wind material might escape to the IGM. In this case, if the 
wind detected in M100 were representative of similar phenomena in other galaxies 
with low to moderate activity, the current estimates of metal and dust content of the 
IGM might be drastically underestimated. 
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1 INTRODUCTION 

Galactic winds are one of the main mechanisms by which 
galaxies expel metals, dust, cosmic rays, magnetic fields, 
etc, to the intergalactic medium. This feedback makes winds 
an important phenomenon in our global understanding of 
galaxies. These winds are usually associated t o galaxies with 
stron g nuclear activity or with starbursts (see lVeilleux et all 
120051 . and references therein) which are powerful enough to 
feed the winds. But very little is known of this phenomenon 
in quiescent, inactive galaxies. In this paper, we explore this 
phenomenon in a relatively quiet galaxy. M100 is a nearby 
(D=16.1 Mpc) galaxy with low nuclear ac tivity (cataloged 
as transition LINER/HII bv lHo et alJll997l ) and a moderate 
star formation rate of a few solar masses per year (accord- 
ing t o its IR luminosity of Lir = 1O 1O,39 L0, Isanders et al.1 
l2003l h Despite this apparent quiescence, we will show that 
this galaxy presents an unexpected fast wind emanating 
from its nuclear region. 

Winds have been studied in highly inclined galax- 
ies where the extraplanar gas can be directly seen 



through imaging of the em ission lines of the ionized 
gas (see IVeilleux et al.l l2005h . This technique has the in- 
convenient of lacking information on the kinematics of 
the outflow. Alternatively, winds have also been studied 
through detection in the spectra of a blueshifted compo- 

ne nt in some interstellar absorption lin es (s ee for exam- 

ple iHeed no recibe p aauctcsckman et al 1 l200d : iRupke et al.l 
120051 : lMartmll200a ). This latter approach lacks spatial res- 
olution, and consequently does not provide information on 
the wind geometry. We have taken a new approach by us- 
ing integral field spectroscopy, which combines spectroscopic 
and imaging techniques. We are therefore able to detect ab- 
sorbing material throughout the field of view and to produce 
a map of the distribution and kinematics of the absorbers. 
In addition, we have been able to discriminate kinemati- 
cally the ionized component of the outflowing material and 
to produce maps for the distribution and kinematics of this 
hotter phase of the wind. This approach provides a unique 
simultaneous global view of these two phases of a galactic 
wind. 
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2 OBSERVATIONS AND DATA ANALYSIS 

The data analyzed in this paper were obtained on 
2002 March 16 at the Observatorio del Roque de los 
Muchac hos on La Palma with the fibre system INTE- 
GRAL l|Arribas et all 119981 ) at the William Herschel Tele- 
scope. The o bservational setup and basic data reduction 
are detailed in ICastillo- Morales et all (|2007l ). We only men- 
tion here the aspects more relevant to the present work. 
The observed spectral range (5600-6850 A) contains a few 
strong interstellar emission lines (Ha, [Nil] AA6548, 6584, 
[SII] AA6716, 6731) as well as several stellar absorption lines 
(NaD doublet and weaker, mostly blended Fel and Cal 
lines) . 

A first inspection of the spectra in the fibres with the 
highest S/N reveals an unusual ratio of the NaD absorp- 
tion doublet with respect to the other stellar lines, with 
th e NaD lines being deeper t han expected (see Fig. 1 
in ICastillo-Morales et al1l2007l ). The most straightforward 
explanation for this is contamination by interstellar NaD 
absorption. We therefore proceed to separate the stellar 
and interstellar components in the fibre spectra. We use a 
synthetic stellar population of 1 Gyr with solar metallic- 
ity to subtract the stellar c omponent of the spectra (see 
ICastillo-Morales et al.l 120071 . for a detailed explanation). 
This synthetic spectrum is used as a template to fit the 
spectrum of each fibre. The free parameters in the fit are 
the velocity, velocity dispersion and (global) line strength. 
It is important to remove from the fitting the wavelength 
range containing the NaD doublet, as it is likely to be con- 
taminated with interstellar absorption. The result of this 
procedure reproduces very well the stellar features in the 
spectra except for the NaD lines. The stellar template fits 
the redmost part of the NaD lines, leaving a strong inter- 
stellar residual which is clearly blue shifted with respect to 
the s tellar component (see Fig. 1 in ICastillo-Morales et al.l 
2007). The presence of blueshifted interstellar absorption 
lines is an u nambiguous signat ure of outflowing material (see 
for example lRupke et al.ll2005l . and references therein) . After 
subtraction of the stellar component, a spectrum of the in- 
terstellar material remains (with both, emission and absorp- 
tion lines). We analyze these interstellar spectra throughout 
the rest of this paper. 



MORPHOLOGY AND KINEMATICS OF 
THE NEUTRAL NA I PHASE 



Most works on the 
ti c winds by means 



study of outflows or galac- 
of the NaD abso rption (e.; 



Heed no recibe paquctcsckman et al 



20001; iRupke et all 



Martini 120051 . 120061 ; 



20021; ISchwartz fc Martini 120041 : 
Rupke et alj l2005h do either select samples of objects 
where the NaD doublet is dominated by the interstellar 
component, or try to make a global rough correction for the 
stellar contribution. We present here a spatially resolved 
separation of these two components. With this aim, we 
apply the analysis described in section [2] to the spectra of 
the 2D collection obtained from the central region of M100, 
deriving full 2D maps of the distribution and kinematics of 
the absorbing material. A first result of this study is that 
regions where the NaD is dominated by the interstellar 




Figure 1. Composite RGB image of the equivalent width of the 
interstellar NaD doublet (blue), K band image (green) and CO 
intensity map (red). Intensities are arbitrarily scaled. The FOV 
of SB3 (outer box corresponding to size of Figs. 2 and 3) and of 
SB2 (inner box corresponding to size of Figs. 6 and 7) have been 
marked. 



component coexist in the same galaxy with regions where 
the stellar component is dominant and there is very little, if 
any, cold interstellar absorbing gas. We would like to point 
out at this point that 2D spectroscopy has been critical 
in detecting the wind thanks to its ability to isolate the 
locations where interstellar absorption is relevant. Indeed, 
the wind might have been unnoticed in a spectrum taken in 
a large aperture where stars dominate the NaD absorption. 

In Fig.[T]we present the equivalent width map for the in- 
terstellar NaD doublet, which ranges between 800 and 2100 
mA throughout the field of view (lower values of the EW 
have been masked out) . The distribution of N al in the map is 
clear ly elongated along the bar position angle l|Knapen et al.l 
1995) with a central maximum and two knots. The cen- 
tral maximum coincides with the nuclear starburst and the 
location of the two knots matches that of the NIR knots 
corresponding to active star forming regions. These knots 
are also related to two strong concentrations of absorbing 
material at these locations that a re clearly seen in the CO 
map (from ISakamoto et al.l Il995i. With the exception of 
these knots, most of the gas responsible for the NaD ab- 
sorption is concentrated in the region interior to the star 
formation ring. The determination of the column density 
of Nal from the spectra is, in general, not straightforward 
(see [Martin 200C], and references therein). In the case under 
study, the situation is worsened because the stellar compo- 
nent must be removed, which can seriously affect the shape 
of the final spectra. Thus, the fi tting of optical dep th and 
covering factor as performed by IRupke et al.l l|2005h might 
not be fully reliable for all spectra (particularly but not ex- 
clusively for those with relatively low S /N ratio) . The equiv- 
alent width is much less affected by the above-mentioned 
effects and is therefore a more reliable estimator of the 
column density. We then estimate the column density of 
Nal from the equivalent width of the NaD doublet using 
a curve of growth technique. We have taken a hybrid ap- 
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Figure 2. Map of hydrogen column density determined from the 
Nal absorption. Contours are spaced every 10 20 cm~ 2 . 
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Figure 3. Velocity field of the absorbing Nal gas. Contours are 
spaced every 20 km/s. 

proac h for this purpose. I n first place, we use a direct fitting 
as in iRupke et alj {2005) on a selected subset of high S/N 
spectra to estimate some average parameters of the redder, 
AA5896, line: covering fraction, Cf = 0.17±0.05, Doppler pa- 
rameter, b = 119±35km/s, and optical depth, to = 0.7±0.3. 
In a second step, we use these average values and the curve 
of growth to find the following relation between the Nal 
column density and the equivalent width of the NaD dou- 
blet: N(Nal) « EW x 2.88 x 10 13 cm~ 2 . Finally, to convert 
N(Nal) to hydrogen column density, we adopt the relation: 
Iog[ N(NaI)l = /31og[N(H )1 +7with/3 = 2 11 an d 7 = -31.3 
from lStokes et all ljl978h as lRupke etall (|2002T ). The result- 
ing map of hydrogen column density is shown in Fig. [2] From 
this map we have calculated the total hydrogen mass in the 
wind to be ~ 2.4 x 10 6 M Q . The velocity at each fibre is 
calculated by fitting two gaussians to the NaD doublet in 
the (stellar removed) spectra. Errors in the velocity deter- 
mination range from 20 km/s to 70 km/s, with the highest 
errors occurring at the locations with the lowest EW. The 



velocity map of the absorbing gas (Fig. [3} presents a rota- 
tion pattern with a major kinematic axis orientation and 
velocity amplitude t hat clearly resembles that of the ionized 
gas in the disk (see ICastillo- Morales et al"1 [2007). However, 
this rotation pattern is globally blueshifted with respect to 
the underlying ionized gas in the galactic disk (v sya = 1567 
km/s) by ~ 115 km/s. The detection of this rotation pat- 
tern in the absorbing gas indicates that this neutral phase 
is very likely disk material that has been dragged upwards 
by the galactic wind, but that still keeps most of its angular 
momentum. It is probably material near the galactic disk, 
on the interface between the fast, hot gas filling the cavity, 
and the underlying quiescent gas in the galactic disk. 

In order to measure the motion of the absorbing gas 
with respect to the underlying galaxy, we have subtracted 
the velocity field of the i onized gas in the disk (see 
ICastillo-Morales et alj 120071 ) from the NaD velocity map. 
The resulting projected velocities range from ~ -210 km/s 
to roughly -55 km/s with an average of -115 km/s. 



4 MORPHOLOGY, KINEMATICS AND 

PHYSICAL CONDITIONS OF THE IONIZED 
GAS PHASE 

The distribution and kinematics of the ionized gas in the 
circumnuclear region of M100 has b een extensively studied 
by se veral authors (see for example Castillo-Mor ales et al.l 
120071 . and references therein) . However, after our finding of a 
strongly blueshifted component in the neutral gas traced by 
the NaD lines, it is reasonable to wonder whether this out- 
flowing gas has an ionized counterpart, and if this is the case, 
whether we can detect it. This detection has a major diffi- 
culty with respect to the neutral gas, namely, the fact that 
we see the outflowing material in front of the bright back- 
ground galaxy. Fortunately, integral field spectroscopic data 
provides kinematical information, and we may try to kine- 
matically di scriminate these two ionized com ponents (see 
for example lArribas fc Mediavillal Il993l . Il994h . 

Visual inspection shows indeed a blue shoulder in the 
emission lines in many spectra taken with the SB3 and SB2 
bundles (see, for instance, Fig. 2} . Unfortunately, the spec- 
tral resolution with the SB3 fibre bundle is not high enough 
(4.8 A), and this blue shoulder is only visible in a few fi- 
bres. The SB2 bundle offers a much better (2.8 A) spectral 
resolution, and the line asymmetry is visible in many more 
fibres. We therefore perform this kinematical discrimination 
in the SB2 bundle only. In order to study this blueward 
component we try to fit two kinematical components to the 
emission lines of each spectra. To do this we simultaneously 
fit two sets of five spectral lines (Ha, [Nil] AA6548, 6584, 
[SII] AA6716, 6731). We assume that, for each set, all these 
five lines have the same velocity and velocity dispersion. The 
theoretical 3:1 line ratio is fixed for the [Nil] doublet. Fig. 
[4] shows the fit to a spectrum where the double component 
is clearly detected. A kinematical double component is said 
to be detected only if the five spectral lines are detected in 
the two sets and the reduced \ 2 of the fit shows a signifi- 
cant improvement of at least 50% with respect to the fit to 
a single kinematical component. As a final check, each fit is 
inspected visually and in case of doubt a single component 
is kept. When this procedure is finished, there is still some 
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Figure 4. Plot of fibre #89 where the detection of a double kine- 
matical component can be seen. Long and short dashed lines show 
the best fit to the disk and wind material respectively. Continous 
line represents the total fit. 
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Figure 5. Composite RGB image of Ho intensity of the wind 
component (red), HST/ACS F555W (green) and Ho intensity 
of the underlying galaxy disk (blue). Intensities are arbitrarily 
scaled. 

ambiguity at some positions on which of the two compo- 
nents belongs to the background galaxy and which to the 
outflow. Continuity in the produced maps (particularly in 
the velocity field of the background galaxy), and the higher 
velocity dispersion of the secondary component are used to 
discriminate between the two. In Figs. [5] |6] we present the 
intensity and velocity maps corresponding to the wind com- 
ponent of the ionized gas phase. The spatial distribution of 
the ionized gas wind component is concentrated around the 
nuclear starburst (although somewhat elongated along the 
major axis) with a bright ridge that extends towards the 
SW and a few weaker patches. The central and SW ridge 
may be enhanced emission from the walls of the expanding 
shell seen in projection. The ionized component of the wind 
fills the region interior to the ring. If there is indeed an ex- 
panding shell in this region it could be possible that the star 
forming ring is somehow related to the wind. 

The velocity field of the ionized gas shows no rotation 
pattern. Nevertheless, the velocity field of the underlying 
ionized gas in the galaxy disk has been subtracted to view 
the motion of the wind with respect to it. It is blueshifted 
in most locations if we make exception of a redshifted re- 
gion in the N. The maximum blueshifted velocity is located 
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Figure 6. Velocity field of the ionized component of the wind 
with respect to the underlying galaxy. Contours show the distri- 
bution of galactic Ho. 

in the nucleus and on a ridge extending from the S to the 
NW, which can be interpreted in terms of the projection of 
the walls of an outflowing cone or shell. These results are 
compatible with a bi-polar outflow (in which case we could 
identify the region of receding velocities with the dimmed 
wind on the far side of the galaxy). 

The wind component has [NII]/Ha ratios in the 1.2 
to 2.7 range and [SII]/Ha in the 0.4 to 1.5 range, typical 
of shock ionization. The galactic component shows lower 
values, with [NII]/Hq of 0.9 for the nucleus, and 1.1 < 
[Nil] /Ha < 1.4 (0.3 < [SII]/Ha < 0.7) for the region in- 
side the star forming ring and of about 0.2 in this ring (see 
Fig. [7}. This indicates that even for the galactic component, 
shock ionization is not negligible in the nucleus, and is even 
more important in the region inside the star forming ring. 
This region also shows the highest velocity dispersion (up 
to 100 km/s) of the region. Whether a hidden wind compo- 
nent sharing the disc kinematics is responsible for the high 
line ratio and velocity dispersion in the galactic component 
in this region around the nucleus, remains unclear. On the 
contrary, the ionization in the stellar formation ring appears 
not to be directly related to the wind. The velocity disper- 
sion of the ionized gas in the wind is in the range between 
150 and 230 km/s, while for the galactic component this 
value is much lower, ranging between 50 and 100 km/s. 



5 DISCUSSION 

In previous sections we have reported the discovery of an 
important wind in the central region of M100. This is cer- 
tainly an unexpected finding, particularly because its extent 
(~ lkpc) and velocity (~ 200km/s) are quite large for a 
galaxy with a moderate nuclear and/or star forming activ- 
ity. In fact, the velocity of this wind is comparable to the 
values found by iMartinl (2006) for ULIRGs. We have es- 
timated the total mass of gas in the wind to be of about 
2.4 x 1O 6 M0. Although this value is abo ut an or d er of mag- 
nitude lower than the values found by iMartfn] (|2006t ). we 
must keep in mind that ULIRGs are forming stars one hun- 
dred times faster. This comparison makes the wind in M100 
even more spectacular and surprising. 
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Figure 7. [NII]/Hq line ratio for the wind (left) and galactic (right) components. Contours show intensity of galactic Ha emission. 



Although the geometry of the wind is uncertain, our 
findings indicate that the galactic wind is most likely driven 
uniquely by the nuclear starburst. The presence of rotation 
in the kinematics of the absorbing neutral gas indicates that 
it is likely disk ISM near the plane that has been entrained 
by the wind at the interface with the hot expanding shell. 
The fact that we have not detected a rotation pattern in 
the wind component of the ionized gas, together with the 
high [Nil] /Ha and [SII]/Hq ratios, indicates that this is 
shocked gas at the top and/or walls of the expanding shell. 
This scenario is also fully compatible with the irregular ge- 
ometry of the ext ended continuum radio emission found by 
iFilho et al.l ||2000| ) , and the diffuse X ray emission reported 



by Immler et al.l (| 1998T ) . This new scenario also opens the 
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possibility that the nuclear star forming ring in M100 is re- 
lated to this wind, either by direct compression of the ISM 
at the edges of the shell, or by mater i al of the wind falling 
down to the disk. In fact, iRvder et all <|200ll ) found evidence 
for sequentially triggered star formation in the ring of M100, 
compatible with this scenario. 

Taking into account that we find projected ve- 
to ~ 210 km/s and velocity dis- 
to 220 km/s for the ionized gas, 
negligible fraction of the gas (with 
* \/2v c (l + ln(r max /r)) 1/2 ~ 500km/s) 
which might be able to escape the gravitational potential 
of the galaxy (take n from a singular isothermal sphere as in 
iRupke et all (|2002h ). Although this rough estimate should 
be taken with caution, it indicates that winds of this type in 
the abundant galaxies with low nuclear activity may be very 
important in polluting the IGM with dust and metals. The 
relevant question then is: how representative is the wind 
in M100 of similar phenomena in other normal, relatively 
unactive galaxies?. This question can only be answered by 
observing a significant sample of low to moderate activity 
galaxies with the techniques used in this paper. 
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